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Abstract 

The aim of this paper is to study a whole class of first order differential inclusions, which 
fit into the framework of perturbed sweeping process by uniformly prox-regular sets. After 
obtaining well-posedness results, we propose a numerical scheme based on a prediction- 
correction algorithm and we prove its convergence. Finally we apply these results to a 
problem coming from the modelling of crowd motion. 
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1 Introduction 



The study of first order differential inclusions started in the 1960s with the theory of the maximal 
monotone operators (see e.g. [5]). Later J. -J. Moreau considered a problem involving a time- 
dependent multivalued operator in [21]. He dealt with the first sweeping process by convex sets 
C(t) included in a Hilbert space: 

^(t)G-a/c(t)(q(t)), (1) 

where die is the subdifferential of the characteristic function of a convex set C. Such a situation 
may be visualized as a point q(t) moving inside C{t) and being pushed by the boundary of this 
convex set when contact is established. This problem can also be written : 

^We-N(C(t),q(t)), 

where N(C, q) is the proximal normal cone to C at q (see Definition 12. ip . He developped a 
so called catching-up algorithm to build discretized solutions and so proved the well-posedness 
of ([I]) under some assumptions on the set-valued map C(-). More precisely, in considering 
some subdivision {Jk)k of the time-interval, the set-valued map C is approached by a piecewise 
constant multifunction taking value Ck on J^. The associated discretized solution q defined by 

yt G Jfc+i , q(t) = qfc+i = Pcfc+i(qfc), 

with qo fixed to the initial value, converges to the unique solution of ([T]). 

Since then, important improvements have been provided by weakening the convexity assumption 
and by considering a perturbed version of this problem: 

^(t)+N(C(t),q(t))9f(t,q(t)). (2) 

In [23], M. Valadier studied sweeping process by complements of convex sets in finite dimension 
without perturbation. Perturbations (even multivalued perturbations) have been taken into ac- 
count by C. Castaing, T.X. Due Ha, M.D.P Monteiro Marques and M. Valadier in [HI [7]. 

Then the main concept, which appeared to weaken the convexity assumption of sets C{t), is 
the notion of uniform prox-regularity. A set C is said to be rj-prox-regular (or uniformly prox- 
regular with constant rj) if the projection onto C is single- valued and continuous at any point 
whose the distance to C is smaller than r]. Under this assumption, the sweeping process without 
perturbation was firstly treated by G. Colombo, V.V. Goncharov in [TOj, by H. Benabdellah 
in [T] and later by L. Thibault in |22j and by G. Colombo, M.D.P. Monteiro Marques in [llj . 
The perturbed problem was later studied by M. Bounkhel, J.-F. Edmond and L. Thibault in 
[3l[221[l3l[Il]. 

In this paper, we consider perturbed sweeping process by a set-valued map Q : [0, T] ^ M'^ 
satisfying that for every t E [0,T], Q{t) is the intersection of complements of smooth convex 
sets. Let us first specify the set- valued map Q. For i £ {1, ... ,p}, let gi : [0, T] x R'^ ^ M be a 
convex function with respect to the second variable. For every t € [0, T], we introduce the sets 
Qi{t) defined by: 

Q,(t) :={qGM^ <7.(i,q) >0}, (3) 
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and the feasible set Q{t) (supposed to be nonempty) is 



Q{t) ■.= f]Qiit). 



(4) 



1=1 



The associated perturbed sweeping process can be expressed as follows: 

^(t) +N(Q(t),q(t)) 3 f(t,q(t)) for a.e. t G [0,T] 
q(0) = qo G Q(0). 



(5) 



This differential inclusion can be thought as follows: the point q(t), submitted to the pertur- 
bation f(t, q(t)), has to stay in the feasible set Q{t). To obtain well-posedness results for @, 
we will make the following assumptions which ensure the uniform prox-regularity of Q{t) for all 
t £ [0,T]. We suppose there exists c > and for all t in [0, T] open sets Ui{t) D Qi{t) verifying 



dH{Qi{t),R''\Ui{t))>c, 



(AO) 



where dn denotes the Hausdorff distance. Moreover we assume there exist constants a, /3, M > 
such that for all t in [0, T], gi{t, •) belongs to C'^{Ui{t)) and satisfies 



VqG C/i(t), a< |Vq5i(t,q)| < A 

VqeM^ |at5i(t,q)| </3, 
VqGC/i(t), |atVq5»(t,q)| <M, 



and 



VqG C/i(t), \Dlg^{t,c^)\ <M. 
For all t G [0, T] and q G Q(t), we denote by I{t, q) the active set at q 

/(t,q) :={ie{l, ...,p}, g,{t,ci)=0} 

and for every p > 0, the following sets: 

Ip{t,(l) :={iG {1, 5i(i,q) <p}- 

In addition we assume there exist 7 > and p > such that for all t G [0, T], 



VqGQ(t), VAi >0, ^ A,|Vq5i(t,q)| <7 

je/p(t,q) 

We will use the following weaker assumption too: 



VqGQ(t), VA, >0, J2 ^^\^^9i{t,fl)\<7 

iel{t,q) 



AiVq5(i(t,q) 

je/p{t,q) 



X] ^iVq5ri(t,q) 

ie/(t,q) 



(Al) 

(A2) 
(A3) 

(A4) 
(6) 

(7) 
(A5) 



(A5' 



In particular, this last assumption implies that for all t, the gradients of the active inequality 
constraints Vqgi(t,q) are positive-lineary independent at all q G Q(t), which is usually called 
the Mangasarian-Fromowitz constraint qualification (MFCQ). Conversely the MFCQ condition 
at a point q yields a local version of Inequality (|A5'p . 
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The notion of uniform prox-regularity allows to adapt the catching-up algorithm because the 
projection onto a uniformly prox-regular set is well-defined in its neighbourhood. However, 
from a numerical point of view, it may be difficult to perform this projection. Here, we study 
a numerical scheme avoiding this difficulty, which is adapted from the one proposed in [15\ \W[ 
118] . The idea is to replace Q{t) with a convex set Q{t,q) (depending on the position). This 
substitution is convenient because classical methods can be employed to compute the projection 
onto a convex set. Yet this replacement raises some difficulties for the numerical analysis which 
are solved in proving that for every q G <5(t), the set Q{t,q) is a good local approximation of 
Q{t) around q. 

The paper is structured as follows: In Section [21 we describe the mathematical framework to 
study the differential inclusion ([5]). By justifying that the set-valued map Q is Lipschitzian 
and takes uniformly prox-regular values (Propositions 12.111 and 12. 9p , we obtain well-posedness 
results for ([5]) in Theorem I2.12[ Then in Section [3] we propose a prediction-correction scheme 
()22p and prove its convergence in Theorem 13.41 Finally we apply these results in two examples 
in Section [H The first situation is a case in point and the second one comes from the modelling 
of crowd motion in emergency evacuation. 

2 Mathematical framework and well-posedness results 
2.1 Preliminaries 

Firstly we recall some definitions and properties to specify the mathematical framework. Here 
we consider a finite dimension space M'^ equipped with its Euclidean structure although these 
notions have been extended in a Hilbertian context. For more details, we refer the reader to 

PS!- 

Definition 2.1. Let S be a closed subset ofW^. 
The proximal normal cone to S at x is defined by: 

N(5, x) := |v G M"*, 3a >0, x G Ps{^ + av)| , 

where 

Psiy) ■= {z G S, ds{y) = |y - z|}, with ds{y) := inf |y - z| 
corresponds to the Euclidean projection onto S. 

Note that for all x G int{S) := S\ dS, N(5, x) = {0}. This concept extends in a certain way the 
notion of normal outward direction for a smooth manifold, as specified in the next proposition. 

Proposition 2.2. Let S be a closed subset o/M'^ whose boundary dS is an oriented hyper- 
surface. For each x G dS, we denote by z^(x) the outward normal to S at x. Then, for each 
x G dS, the proximal normal cone to S at x is generated by ^{x), i.e. 

N(5,x) =M+i/(x). 



We refer the reader to Proposition 3.2 of [24j for a detailed proof. 
Using [9], we recall the concept of uniform prox-regularity as follows: 
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Figure 1: Ty-pr ox-regular set. 



Definition 2.3. Let S he a closed subset ofW^. S is said rj-prox-regular if for all x G dS and 
V G N(5', x), |v| = 1 we have: 

B{x + riv,ri) n 5 = 0. 
Equivalently, S is r]-prox-regular if for all y £ S, x £ dS and v G N(S', x), 

(v,y - x) < |^|y - xp. 



In other words, S is 77-prox-regular if any external ball with radius smaller than r] can be rolled 
around it (see Fig[T]). Moreover, this definition ensures that the projection onto such a set is 
well-defined and is continuous in its neighbourhood. Note that a closed convex set C C M"^ is 
00-prox-regular. We claim a technical lemma which will be useful later. 

Lemma 2.4. Let S be a closed convex set included in W^. Then for all x £ S and w G M'^.' 

w G N(5, x) ^ x = Ps(x + w) (a) 

^ Vy gS, (w,y-x) <0 (6) 

^ G M'^, (w,^) < |w| (i5(^ + x) (c) 

3k>0, 36> 0, Vv G M'^, |v| < 5, (w, v) < A; dsiv + x) (d) 



The elementary proof is left to the reader. 



2.2 Uniform pr ox-regularity of sets Q(t) 

This subsection is devoted to specify the proximal normal cone to Q{t) (defined by ^) for 
t £ [0,T] and to justify the uniform prox-regularity of this set under the assumptions (|Aip . (|A4p 
and (1A5'|) . We point out that the time-dependence plays no role in this study. 
First we deal with only one constraint. Fix i G {1, . . . ,p} and consider the set Qi{t) defined by 
([3]). The smoothness of the function gi{t,-) allows us to apply Proposition 12.21 and to deduce 
the expression of the proximal normal cone to Qi{t). 

Corollary 2.5. For all t G [0,r] and q G dQi{t), 

N(Q,(t),q) = -M+Vq5i(i,q)- 
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By Definition 12. 31 the constant of prox-regularity equals to the largest radius of a "rolling external 
ball". In order to estimate its radius, tools of differential geometry can be used. More precisely, 
to show that the set Qi{t) is uniformly prox-regular, we can apply the following theorem proved 
in [EJ. 

Theorem 2.6. Let C be a closed convex subset o/M'^ such that dC is an oriented hypersurface 
ofM.'^. We denote by I'ci'^) the outward normal to C at 'k and by pi(x), .., (x) > the 
principal curvatures of C at x. We suppose that 

p := sup sup Pi(x) < oo. 
xg9C l<i<d-l 

Then S = M.'^ \ int{C) is a rj -prox-regular set with rj = —. 

Proposition 2.7. For all t G [0, T] and i G {1, . . . Qi{t) is r]Q-prox-regular, with r]Q = 

Proof. Let t G [0,T], by Assumption ()A1|1 . the complement of Qi{t) is obviously the interior 
of the closed convex set C = {gi{t, •) < 0} which satisfies the assumptions of Theorem 12.61 
The constant of prox-regularity of Qi{t) can be obtained by calculating its principal curvatures 
which are the eigenvalues of the Weingarten endomorphism. Let q G dQi{t) = dC, the outward 
normal to C at q is equal to — i/(q), where 



|Vq5i(i,q)r 

We can specify the expression of the Weingarten endomorphism Wq. For every tangent vector 
h G T^{dQi{t)), 

Wq(h) := -DKq)[h] = T^r^J—T^Pi.i^))^ (D25,(t,q)[h]) , 

with 

P;,x(x) = X - (x,i/)z^. 

By ([XT]) and (jUl) . for ah q G dQi{t) C Ui{t), the eigenvalues of Wq are bounded by Ma^^, 
which ends the proof according to Theorem 12.61 . □ 

Now let us study the feasible set Q{t) that is the intersection of all sets Qi{t). We begin to 
determine its proximal normal cone. 

Proposition 2.8. For all t G [0, T] and q G Q{t), 

N(Q(t),q) =^N(Qi(t),q) =- K+Vq5^(i,q). 

«G/(t,q) 



We refer the reader to Proposition 2.16 in [T7j for a detailed proof. We emphasize that the 
inclusion J2^iQi(t),<l) C N(Q(t),q) always holds. However, the given proof of the other in- 
clusion requires the convexity of functions gi{t, •) and the boundedness of gradients Vqgi{t, •) 
(Assumption (jATil ). 

We now come to the main result of this subsection: the set-valued map Q takes uniformly prox- 
regular values (with a time-independent constant) which rests on the inverse triangle inequality 
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Proposition 2.9. For every t G [0, T], Q{t) is rj-prox-regular with 

rjo a 



rj 



7 M-f' 



Proof. Let t E [0, T] and consider q G dQ{t) and v G N(Q(t), q) \ {0}. By Proposition ESI there 
exist Oi >0 such that 

v = - ^ aiVqgi(t,q). 

ie/{t,q) 

By definition Q of /(t, q), q G dQi{t) for every i G /(t, q). By Proposition [221 for all i G /(t, q), 
— OiVq q) G N(Qi(t),q), so we have by Definition! 



(-aiVqfifi(t,q),q- q) < 



2% 



q - qp, Vq G 



Since Q{t) C Qi{t), by summing these inequalities for i G /(t,q), we obtain 



^ aiVqC/i(t,q),q- q\ < ^ 

je/(t,q) / «e/(t,q) 



75 q - q . Vq G Q{t). 

2r/o 



^Prom ()A5'P . it follows that 



X] aiVq5(i(t,q),q - q \ < 



7 

2% 



ie/(t,q) 



|q-q|2, VqGQ(t). 



□ 



\ iG/{t,q) / 

We deduce from Definition 12.31 that Q{t) is ry-prox-regular with rj = 
2.3 Lipschitz regularity of Q 

Here we check that t — > (5(t) varies in a Lipschitz way. Aiming that, we need this technical 
lemma: 

Lemma 2.10. There exists 6 > such that for all t G [0,T] and q G Q{t): 
3uGM^ |u| = 1, ViG/p(t,q), (Vq5i(t,q),u) >5, 



(8) 



Proof. We set the following cone 

Np(Q(t),q) :=- J] M+Vq5.(t,q) 

«G/p{t,q) 

and its polar cone 

Cp(Q(t),q) :=Np(Q(t),q)° := {w G Vv G N^(Q(t), q) , (v,w) <0}. 

According to the classical orthogonal decomposition of a Hilbert space as the sum of mutually 
polar cones (see [20]), we have: 



Id = P 



Np(Q(t),q) + PCp(Q(t),q), 
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where P denotes the Euchdean projection. 

So let us consider for i £ Ip{t, q) the corresponding decomposition of V q^gi{t^ q): 

Vq<7i(t,q) =«i+&i eNp(Q(f),q)+Cp(g(t),q). 

Assumption (lAip gives us: \bi\ < |Vq(7i(t,q)| < (3. Since ai G Np((5(t),q), it can be written: 
ai = -Yl '^i^q djit, q) , > involving 



|Vq5j(t,q) - Oil 



Aj' Vq 5i (t, q) + (1 + Ai) Vq giit, q) 



Then using the inverse triangle inequality (jASp and Assumption (lAlh . we get : 



As a consequence, it comes: 



|6^l>-(5]A, + l) >- 



7 



7 



7 



Since £ Np{Q{t),q) and = PNp(Q(t),q)(Vq5fi(t, q)), we obtain 

2(6i,-Vq(7*(t,q)) = |6i - Vqgi(t,q)|2 - \bi\^ - |Vq5i(t,q)P 
= |aiP - l^iP - |Vq5i(t,q)p 
,, ,9 

<-6^'<- — . 

7^ 



Now we set 



I Y.ieip{t,ci) 

This is well-defined because (llOp and Assumption (lAip imply that for any j G Ip(t,q) 



(9) 



(10) 



(11) 



«6/p(t,q) 



E ^-Vq*(t,q))>^(&„Vq<7,(t,q))>^. 



Ue/p(i,q) 

Then (El) follows from & and ([12]) with 



(12) 



a 



272^/3 



□ 

Proposition 2.11. T/ie set-valued map Q is Lipschitz continuous with respect to the HausdorjJ 
distance. More precisely there exists Kl > such that 

Vt,sG [0,T], dH{Q{t),Q{s))<KL\t-s\. 

Proof. Consider t,s £ [0, T] and q G Q{t), let us construct a point closed to q belonging to 
Q{s). Let u given by Lemma 12.101 we introduce z(/i) := q + hu with > 0. We claim that for 
h < hi := min(c, p/(/3 + (5)) , 



Vi G {1, . . . gi{t,z{h)) > h5 
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(where c and 6 are introduced in ()A0|) and ()2.10p ). Indeed for h < c, due to the convexity of 
gi{t,-), it comes 

gi{t,(i+hu) > 5i(t,q) + /i(Vq5i(t,q), u). 
As a consequence, for i G /p(t, q), 

gi{t,c[ + hu) > h6, 
by ([H]). Moreover for every i ^ /p(t,q), according to Assumption ()Aip . 

gi{t, q + /lu) > p — hf3 > h5 

if h < Thus for /i < hi, we have q + hu.) > h6 for ah i E {1, . . . ,p}. That is why we 

phi 

deduce from Assumption (|A2p that z(/i) G (5(s) if /i^ > (3\t — s\. Setting £ := — — , if |t — s| < £, 



it can be written that 

l3 

dQ(s){^} < ^ inf , |q - zWI < -jlt - s\. 

^ ' h5>l3\t-s\ 

Consequently we obtain if |t — s| < i, 



dH{Q{t),Q{s)) = max! sup (iQ(^)(q), sup (iQ(t)(q)) <^|t-s| 



This inequahty is actuahy satisfied for any t,s £ [0,T]. To check it, it suffices to divide the 
corresponding interval into subintervals of length i and to apply the triangle inequality. So the 
required result is obtained with Kl = 

□ 

2.4 Well-posedness results 

We now come to the main result. 

Theorem 2.12. Let T > and f : [0,T] x M"^ ^ M'^ be a measurable map satisfying: 

Vq G Q{s) , f(-,q) is Riemann-integrable on [0, T] (13) 

se[o,T] 

3iC>0,VqG IJ Q(s), VtG [0,r], |f(t,q)-f(t,q)| <i^|q-q| (14) 

se[o,T] 

3L>0,Vqe U Q(s), VtG [0,r], |f(t,q)| <L(l + |q|). (15) 

se[o,T] 

Then, under Assumption IIA5\) for all qo G Q{0), the following problem 

^(t)+N(Q(t),q(t))9f(t,q(t))/ora.e. t G [0,T] ^^^^ 

q(o) = qo, 

has one and only one absolutely continuous solution q satisfying q(f) G Q(t) /or every t £ [0,T]. 

Proof. As the set-valued map Q varies in a Lipschiz way by Proposition [2TT] and for all t G [0, T], 
Q{t) is 7/-prox-regular by Proposition 12.91 we can apply Theorem 1 of |I3]. □ 

We now consider a constant set-valued map Q (the constraints gi are supposed to be time- 
independent). In this particular case, we only require Assumption ()A5'P instead of Assumption 
[5]). 
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Theorem 2.13. Under Assumptions fA 0\) . liA liA4\ ), ^A5]) with time-independent set Q, let 
T > and f : [0, T] x Q ^ M*^ be a measurable map satisfying: 



Vq G Q , f(-,q) is Riemann-integrable on [0, T] 

3if > , Vq G Q , Vt G [0, T] , \{{t, q) - f (t, q)| < K|q - q| 

3L>0, VqGQ, VtG [0,r], |f(t,q)| < L(l + |q|). 



(17) 
(18) 
(19) 



Then, for all qo G Q, the following problem 



^(t)+N(Q,q(0)9f(t,q(i))M a.e. t G [0, T] 
q(o) = qo, 



(20) 



has one and only one absolutely continuous solution q taking values in Q. Moreover the map q 
is also solution of the following differential equation: 



( da 

^(i) + PN(Q,q(t))(f(i,q(t))) = f(t,q(i))M a.e. t e [0,T] 
. q(0) = qo, 

For this last point we refer the reader to Proposition 3.3 in [2]. 



(21) 



Remark 2.14. The assumptions about the perturbation f can be weakened (see I13f). However 
they are required in the numerical analysis of the scheme proposed in the following section. More 
precisely, we need Assumption I^13\) to prove Lemma \3.(k 



3 A numerical scheme 
3.1 Presentation 

We present in this section a numerical scheme to approximate the solution of ()16p on the time 
interval [0,T]. This scheme is adapted from the one proposed by B. Maury for granular media 
in [15]. Let n G N*, h = T/n be the time step and = kh be the computational times. We 
denote by q^ the approximation of q(t^) with qg = qo- The next configuration is computed as 
follows: 

q^+i:=PQ(,.^^,q.)(q^ + M(tLqD) (22) 

with 

Q(t,q) :={qGM^ ffi(t, q) + (Vq q), q " q) > Vf} for q G ^/(t) :=f|^i(*)- 

We recall that all the gradients Vq5(j(t,q) are well-defined provided that q G U{t). The set 
(5(t, q) can be seen as an inner convex approximation of Q{t) with respect to q. This scheme 
is a prediction-correction algorithm: predicted position vector q^ -|- /if(t^,q^), that may not be 
admissible, is projected onto the approximate set of feasible configurations. 
Let us check that this scheme is well-defined for h < with c and respectively given by 
Assumption ()A0|1 and Proposition ()2.1ip : 

Proposition 3.1. Assume that HKl < c, then for all k < n — 1 

qfc+i G Q(ifc+i,qfc) c g(t^i) c u{ti^,). 

Thus every computed configuration is feasible . 
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Figure 2: Theoretical and numerical projections. 



Proof. Since q[f = qo € Q{0) and dniQiO), Q{ti)) < KlH, we get 



ci^ eQ{t^) + hKLBiO,l)cU{t^), 

according to (|AOp and Proposition 12.111 So (5(ti,qo) is well-defined and is included in Q{ti) 
due to the c onvexity of the functions Then as q" e (3(ii,qo), q? e Qit^) C (5(i2) + 

hKLB{0, 1) C U{t2). By iterating, we end the proof. □ 



By Lemma 12.41 (a), it can be checked that 



h 



N(Q(t^+l,q^),q^+l)9f(tLq^) 



(23) 



so that the scheme can also be seen as a semi- implicit discretization of (I16p . where the cone 
^iQ{'tk+i^'ik)''ik+i) stands for an approximation of N(Q(t), q(t)). 



In Figure [21 we illustrate the set Q{t^j^i) C M'^, intersection of sets whose boundaries 

are plotted in solid line. The set Q{t^^i, q^) is delimited by the dashed line. The theoretical and 
numerical projections, respectively q^^;^ := PQ{t^ i)(^k ~^ ^^(^^5^^)) s-i^d q^^;^ are represented 
(for two examples of f (t^, q^))- Indeed, as Qif^^i) is uniformly prox-regular, the projection onto 
Q{t^+i) of q^ + /if(t^,q^) is well-defined for h small enough. The replacement of (5(ifc+i) by 
the convex set Qitk^i,fl^) is convenient because it allows us to use classical numerical methods 
to compute this projection. However, it raises some difficulties to prove the convergence of the 
scheme. On the one hand, we have to check that this approximation is sufficiently accurate (see 
Lemma 13. 2p . On the other hand, the set Q{t,q) does not vary smoothly enough to directly 
apply results about sweeping process (this point will be specified in Remark l3.7p . 
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Lemma 3.2. For all t G [0,r] and q G Q{t), 

N(Q(t),q) =N(Q(t,q),q). 



Proof. We need the well-known following result (see e.g. [8]) which rests on the reformulation 
of the minimization problem under constraints (associated to the projection), in terms of a 
saddle-point problem. 

Lemma 3.3. Let t G [0,T], q G Q{t) and q G M'^. Then, p = q)(q) "is equivalent to the 
existence of \ G (M"^)p such that (p, A) satisfies: 

yi, gi(t,q) + (Vq(7i(t,q),p-q) >0, (Pt,,,q-) 

Yj^i (9i(*,q) + (Vq<7i(i,q),p-q» = 0. 



First we check that N(Q(t),q) C N(Q(t,q),q). 

Let w G N(Q(t),q). According to Proposition 12.81 this vector can be written 

w = - l^i'^n9i{t,fl) 

iel{t,q) 

with nonnegative reals Hi. So by definition of Q{t, q), we get for all p G Q{t, q), (w, p — q) < 0. 
Applying Lemma 12.41 (b). we deduce that w G N((5(t, q), q). 
It remains to prove N((5(t, q),q) C N{Q{t),q). 

Letting w G N((5(t, q), q), we can write by Lemma 12.41 (a), q = -Pq^^ q) (q + w) which implies by 
Lemma 13. 3[ 

q= (q + w) + ^^iVq5'i(i,q), where (q,/i) G Q{t) x (M+)p satisfies System (Pt^g^g+w)- 

Consequently, w = - ^ gi{t, q) , /Xj > 0. The third relation in (Pf^g^g+w) is Yl f^idiit, q) = 
0, so Hi = if 5i(t, q) > 0, which concludes the proof. □ 

3.2 Convergence 

Before stating the result of convergence, we introduce some notations. We define the piecewise 
constant function f" as follows, 

r{t) = f(t^,q^) if t G [tl,tl+il k < n and r(T) = f(t^i,q^i). 

We denote by q" the continuous, piecewise linear function satisfying for G {0, . . . , n}, q"(t^) = 
q^. To finish, we introduce the functions p and 9 defined by 

p-{t) = tl and r (t) = tl^, if t G [tlXk+il p'^iT) = T and ^ (T) = T. 

Theorem 3.4. With the assumptions of Theorem \2.1^ q" tends to q in C^([0, T], M'^), where 
1 1-^ q(t) is the unique solution of ( figj) . 

Proof. First we prove the boundedness of the sequence (q") and to do this, we use the same 
arguments as developed in p3] with the following lemma. 
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Lemma 3.5. There exists D > such that for h small enough, 

'^Q{q+„<iV^'ik) < Dh. 



Without loss of generality, we can assume that T < 1/{2D + 4L) and h small enough to apply 
Lemma 13.51 

It follows from (1151) that 



<h{D + L){l + \cil\). (24) 



We deduce that 

|q^i-q^|<Ml? + 2L)(l + |q^|) (25) 

which implies that 

|q^ll<|q^^|+Ml? + 2L)(l + |q^|). 

Hence 

n-l 

\cC\ < |qo| + h{D + 2L) V(l + |q^|) < |qo| + nh{D + 2L)(1 + max |q^|). 



0<fc<n 

A:=0 ~ ~ 



Since T = nh, we deduce that 



As a consequence, we obtain by letting C := 2(|qo| + T{D + 2L)) 

l|q"l|oo < C. (26) 

By (gSl) and (IMl), we have 

^ <{D + 2L)il + C). (27) 



dt 

By applying Arzela-Ascoli's Theorem, it can be easily shown that there exists a function q G 
W^'°^{[0,T],Q) and a subsequence (still denoted by q*^) satisfying 

q" > q uniformly in [0, T]. 

n^oo 

By (f26]l . it can be specified that 

llqlloo < C. (28) 

Since the time-interval is bounded, we get 

Now let us check that the limit function q satisfies the differential inclusion ()16p . The beginning 
of the proof uses classical tools developped e.g. in [14J. We want to show that 

^(t)-f(t,q(t))G-N(Q(t),q(0) for a.e. t € [0,T] 
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which is equivalent to 



by Lemma El 
Lemma 3.6. 



^(t) - f {t, q(t)) e -N(Q(t, q(t)), q(t)) for a.e. t E [0, T] 



r-f(-,q(-)) tn L\[0,T],R''). 



This lemma (later proved) implies that 



dt 



dt 



f(-,q(-))in Li([0,r], 



Consequently, by Mazur's Lemma, there exists a sequence G L ([0, T],]R ) satisfying 



and 



E Conv 



dq 



fiq'' 



P, > n 



(29) 



.^-f(,q(.))in L\%T], 



Extracting a subsequence, we may suppose that 

dq 



z 



f(-,q(-)) a.e. in [0,r]. 



Furthermore, Inclusion (1230 can be rewritten for almost every t G [0,r], 
dq" 



dt 



■(t)-r(t)e-N(Q(r(t),q"(p«(t))),q"(r(t))). 



(30) 



(31) 



Let t G [0, T] such that z„(t) tends to z(t) and the above differential inclusion holds. By (p9 
it yields 

G (z"(t),0 < sup /^(t) - f^(t), ^ 

k>n \ OX 



Passing to the limit, we obtain 



G M°', (z(t),^) < limsup 



dq" 



(32) 



From (I3ip and from Lemma 12.41 (c) , we get 



Vn, VI GM^ (^(t)-r(t), I) < 



dq" 



dt 



(t)-r(t) 



d, 



'(3{e"{t),q"{p"{t))) 



:q"(r(t))-i). 



Hence, by (p^ and ([26]) . this inequality can be specified 

Vn, VI G (|^(t) - r (t), 1^ < (Z) + L)(l + C) (q"(r (t)) - |). (33) 

Proposition 13.81 with q = q(t), q^ = ^{p^{t)) and t„ = 0"(t) implies there exists v > Q such 
that for ah | G M"', ||| < v, 



Q(0"(t),q"(p"(t))) 



:qW-0-^^Q(t,q(i))(qW-0 



0. 



(34) 
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Finally, by passing to the limit in (I33p . it comes from (I32p : 

'^'^Q(t,q(t)) 



V^GM^ 1^1 {z{t),i) <{D + L){l + C)do,^,a)Mt)-$), 



which is equivalent to 

z(t)G-N(Q(t,q(t)),q(t)), (35) 



by Lemma 12.41 (d) . The required result follows from (I30p provided we prove Lemmas 13.5 
and Proposition 13. 8[ □ 

Remark 3.7. To estimate the upper limit in i32^) . we give a different method than the one 
written in Jj^[ , \13^ . However the main argument is the same. Indeed in \13^ . moving sets C{t) 
(only depending on time) are considered and the assumed regularity of the set-valued map C 
involves 

dc{t^){v) dc(t){y),yy, when t" t 
which looks like {^)- This property allows the authors to prove Proposition 2.1 of JT^], which 



implies 13^) from ^Wj). In our case, we just need (and only have) a local continuity. More 
precisely, letting 

Q:={(t,q), tG [0,r], qGC/(t)} 

and 

Oe:={(t,q), te [0,T], qGQ(t)}, 

the map 

: S7 X 1 — > M 

(t,q,q) ^Q(t,q)(q) 

is continuous on a neighbourhood of the set {(t,q, q), q = q} in x . In fact, the next 
proposition ensures that the map 

(t,q,q) PQ(t,q)(q) 
is continuous on the same neighbourhood (which is sufficient). 

3.3 Proof of the claimed results 

Proof of Lemma \3.5[ This proof is similar to the one developed for Proposition 12.111 Let u 
given by Lemma [2. lUI with t = and q = q^ G Qi^k)' introduce z(s) := q^ + su with s > 0. 

We claim that for s > (where c is introduced in (lAOh ). 

Vi G Ip(tLqD, A, := gM+i,^l) + {V^gM+i,crk) , zis) - q^) > 0. 

Indeed, by Assumption ()A2p and Proposition 13.11 fl'i(ifc+i,qfc) > giitk,^^) — Ph > —f3h. As in 
Proposition 13.11 it can be proved that for all s G £ Then it follows from 

Assumption ([A3]) that 

(Vq<7,:(i^i,q^), z(.) - q^) > (Vq 5^(t^, qD , z{s) - o^l) - Mhs. 
For h < ^ and by definition of u, 

\ > -Ph - Mhs + sS>s--h(3>0. 
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Moreover, for i ^ /p(t^,q^), Assumption (IA2p yields 

g.{tl^„ce,)>g,{q,cil)-Ph>p-ph. 
Consequently for i ^ Ip{f^,q^) and /i < 

g^{tl^„cll) + (Vq5.(i^i, q^) , z(.) - q-) > p I - 

For /i small enough, for s G [^^, z(s) belongs to q^). So we have proved that 

2/3, 



□ 



Proof of Lemma \3.(A We want to show that f" weakly converges to f(-,q(-)) in L^([0, T], M'^). 
This property rests on the boundedness of the sequence (f")n in -L°°([0, T], M'^). Let us consider 
* G L~([0, r],M'^) supposed to be continuous. 



{r{t),-^{t))dt 

n—\ ^fn n — 1 f.fn 



/ {m,ci{tt)),^{t))dt + Y^ / (f(t^,q"(t^))-f(t^,q(t^)),^(t))dt:=/^+/2. 
|/2|<Kr||M'|U||q"-q|| 



fc=0 ^fc 

By (HI), 



0. 



Furthermore 



n-l 



n-1 



5]Mf(t^,qO,vI/(t^))+j; / '^\f(t^q(t^)),vl/(t)-vl/(t«))dt:= + 

li'^ 

The map f(-,q(-)) is bounded by L(l + C) according to ([15]) and ([28j) . Let us check that 



J, 



n— ++CO 



(f(t,q(t)),v&(i))dt. 



(36) 







For all 7] > 0, there exists m E N such that for all t,s £ [0,T], |q(t) — q(s)| < if |t — s| < 
5 := T/m since q G Tyi'~([0, T], Q). For all t,s G [0,r] satisfying \t - s\ < 6, yields 
|f(t,q(t)) - f(t,q(s))| < Krj. For n > m. 



n-l 



»t-l 



4 = 5] (t^, q(t^)) - m,^{x{tm, ^m) + E Hm,ci{x{tm, ^m), 



where x(0 



fc=o 

: 6 



k=0 



Obviously, 
1 

h{f{tl, q(t^)) - i{tl, q(x(t^))), 



n-l 



k=0 



Moreover 



n — 1 



m—1 



fc=0 



j=0 k 

t^S[i5,(i+l)5[ 
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Since for all i £ {0, ...,m — 1}, f(-,q(i(5)) and ^ are bounded and Riemann integrable by (1130 . 
we have 



'—^ n- 

k 

Thus, as previously, 



iS 



(f(t,q(i<5)),M/(t))dt. 



lim sup 



n— 1 



k=0 



(f(t,q(t)),M/(t))dt 



< TKr]\\^\ 



This concludes the proof of (|36|) . Moreover 

|J^| <rL(l + C7) sup 

\t-s\<T/n 

Then, we deduce that for all continuous map ^ £ I/°°([0, T], II 



n—>+oo 



0. 



n— >+oo 



(f(t,q(t)),M'(t))dt. 



(37) 



Now let <I> G L°°([0, T],M ), there exists a sequence of continuous maps {^p)p satisfying 



0. 



p— ++00 

Let e > 0, there exists p G N such that L(l + C)||^p — <l>||i < e. For that integer p, we write 



< 



< 



{r{t),^t))dt 

T 



(f(t,q(t)),<I>(t))(it 



f {r{t)-f{t,ci{t)),^p{t))dt+ [ (f(t,q(0),^p(t)-f(t))dt+ / {r{t),m-'i'p(.t))dt 

^0 JO 

(p(t)-f(t,q(t)),^p(t))idt+ rmt,<im%it)-m)\dt+ r\nt),m-%m 

Jo Jo 



<Ki + K2 + Ks. 

For i = 2 and 3, iTj < L{1 + C)||^'p — <I>||i < e and by (1370 . we know that for n large enough, 
iCi < e. So we deduce that for all $ G L°°([0,r],] 



{rit),m)dt 



(f(t,q(t)),$(t))dt. 



n— >+oo 

which concludes the proof. □ 

Proposition 3.8. Let t G [0, T], q G Q{t), a sequence (t„) G [0,T]^ converging to t and a 
sequence (q„) G {U{tn))n converging to q such that dg^^ ^ )(ctn) ^e?^rfs to zero. For all q G M*^, 
we denote p (respectively p„j i/ie projection o/q onto (5(t, q) (respectively onto Q{tn, qn) J- T/ien 
i/iere exists u > Q so that for all q G i?(q, i^), i/ie sequence (pn) converges to p. 



Proof. By Lemma [3. 3 1 the vector p satisfies System (Pt,q,q). We obtain similar systems denoted 
by {Ptn,q„,q) for all p„ in substituting t, q and Aj by t„, q„ and A". The following lemma (which 
will be later proved) claims that the nonzero Kuhn- Tucker multipliers A^, A" are associated to 
an index i belonging to /(t, q) (defined by 
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Lemma 3.9. There exist > and Mq G N such that for all n > Mq and all q G S(q, v), we 
have : 

\Vn - q| < 2l/, 

and 

Ar = Ai = 0, i/gi(t,q) >0, 

where (p,A) and (p„,A") are respectively solutions of [Pt^q^q) and iPt„,q„,q)- 

There exists Mi G N such that for n > Mi, |q — q^l < c. For all s G [0,T] and q, x G we 
denote by Ai,^q{s,x) the d x |I(t,q)| matrix defined as follows: 

Let Mq and be fixed by Lemma 13.91 With the previous notation, for n > max(Mo,Mi) and 
q G -B(q, i^), the first equation of {Pt„.q,^,q) can be written as 

^t,q(tn, qn)[A"] = p~„ - q. 

By dJU) and ([AF]) . we have 

|A"|<^IAq(t,q)[A"]|, 
a 

where represents the Euclidean norm of /X G mI^(*''i)I. Moreover 

IAq(i>q)[A"]| < |Pn-q| + |(Aq(in,qn)-Aq(*>q))[A"]|. 

Thus by Lemma |3.9^ 

|A-| <l(2i/ + e„|A"|), 
a 

with £n ■■= ||(^t,q(t„,q„) - At^q{t,q}\\. Hence 



1 _ £!i2^ < ^ 



Furthermore, < -v/P™aXjgj(t q) |Vq(7j(t„,q„) — Vq (7j(t,q)|. Since q„, G it can be written 

|Vqffj(tn,qn) - Vq5rj(t,q)| < | Vq q„) - Vq5rj(t,qn)| + |Vq3j(t,qn) - Vq5(j(t,q)| 

< M(|t-t„| + |q-q„|), 

by Assumptions (jASP and ()A4p . Finally 

en < V^^d* - *n| + |q - qnl) > 

and so we deduce that the sequence (A") is bounded. A convergent subsequence can be also 
extracted and by passing to the limit in System {Pt„,q„,q), we obtain that the corresponding 
subsequence of (p„) converges to a point poo which satisfies System {Pt,q,q)- Lemma [3^ implies 
that Poo = p and we can conclude by compacity arguments. □ 



It remains to prove Lemma 13.91 

Proof of Lemma \3. (A This result is a consequence of the third equation of problems {Pt,q,q) and 
{Pt„,qn,q)- By definition of /(t,q), 

3e > 0, Vi ^ /(t,q), gi(t,q) > 2e. 
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Setting u = as (qn)n and (tn)n respectively converge to q and t, it yields 

3Mo > 0,Vn > Mo, gi{tn,<ln) >e, Vz ^ /(i,q) , '^(3(t„,q„)(qn) < ^ and |q„ - q| < ^. 
Let q E i?(q, z^) , since p = Pg^-^ q)(q) and q G q) (because q G Qit)), we obtain 

|p — q| < |q — q| < 

and consequently, 

|p — q| < |p — q| + |q — q| < 2|q — q| < 2u. (38) 
Moreover, as p~„ = P^^^^ q^-)(q),we get 

|Pn - q| < |qn - q| + C^Q(t„,q„)(qn) < |q " q| + |q " qnl + C^Q(t„,q„)(qn) < Z/ + ^ + ^ = 2l/. 

Hence, 

|Pn - qnl < |Pn " q| + |q " qn| < 2|q - q„| < 4l/. (39) 

For i ^ /(t,q), it follows from ([MD and ([MD 

5i(t,q) + (Vq5i(t,q),p-q) > 2e-2i^/5 > and fifi(t„, qn) + (Vq 5i(t„, q„), p~„-q„) > e-AuP > 0. 

The third equation of problems {Pt,q,q) , {Ptn,qn,q) and the nonnegativity of Kuhn- Tucker mul- 
tipliers Aj, permit us to conclude. □ 

4 Applications 

4.1 A case in point 

Here we deal with a simple case associated to a well-known game: the labyrinth tabletop game. 
The aim of this game is to maneuver a steel ball through a wooden labyrinth by tilting the 
surface in order to move the ball to the target. The angle of the labyrinth need to be carefully 
controlled so that the ball will not fall into the holes. 

We consider a particular labyrinth represented in Figure O with several holes (corresponding to 
the black disks). Initially the ball is located at the starting point S and it has to be rolled to 
the target T. The quasi-static evolution of this ball can be described by the following first order 
differential inclusion: 

^(t)+N(Q,q(t))9f(t), 

where q G is the position of the ball in the labyrinth surface, Q is the complementary set 
of the grey obstacles and the holes and f(t) is the perturbation caused by the leaning of the 
labyrinth plane. In Figure [3l we plotted the trajectory obtained with the following field f: be- 
ginning with f = fi, the ball arrives at the point Pi, then with f = £2, the ball goes to the point 
P2 and finally submitted to f = £3 it reaches the target. 

Obviously, the set Q is uniformly prox-regular and Assumptions (lAlh . (|A4p and (IA5'|) are 
satisfied. The evolution problem is well-posed by Theorem 12.121 and the trajectory of the ball 
can be estimated by the numerical scheme presented in Section O 

4.2 An application to a crowd motion model 

The aim of this subsection is to apply the previous results to a model of crowd motion in 
emergency evacuation. We refer the reader to [HI [T7] for a complete and detailed 
description of this model. 



19 



Figure 3: Labyrinth game. 



4.2.1 Presentation of the model 

We quickly recall the model. Its asset is to directly handle numerous contacts which are char- 
acteristic of emergency situations. This model allows us to deal with local interactions between 
people in order to describe the whole dynamics of the pedestrian traffic. This microscopic model 
for crowd motion rests on two principles. On the one hand, each individual has a spontaneous 
velocity that he would like to have in the absence of other people. On the other hand, the 
actual velocity must take into account congestion. Those two principles lead to define the actual 
velocity field as the projection of the spontaneous velocity over the set of admissible velocities 
(regarding the non-overlapping constraints). 

We consider persons identified to rigid disks. For convenience, the disks are supposed here 
to have the same radius r. The center of the i-th disk is denoted by G M^. Since overlapping 
is forbidden, the vector of positions q = (qi, ..,qAr) G has to belong to the "set of feasible 
configurations", defined by 

Q := {q G M2^, A,(q) > Vi / j} , (40) 

where Djj(q) = |qj — qj| — 2r is the signed distance between disks i and j. 

We denote by U(q) = (f7i(qi), .., f^Af(qAf)) G the global spontaneous velocity of the crowd. 
We introduce the "set of feasible velocities" defined by: 

Cq = {v G R^^, Vi < j Ai(q) = ^ (G,, (q), v) > O} , 

with 

Gi,(q) = VA,(q) = (0,...,0,-ey(q),0,...,0,ey(q),0,...,0) GM^^ 
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and ejj(q) = . The actual velocity field is defined as the feasible field which is the closest 

to U in the least square sense, which writes 

^ = Pc,(U(q)), (41) 

where denotes the Euclidean projection onto the closed convex cone Cq. We deduce from 
Farkas Lemma (see e.g. [8j) and Proposition 12.81 the following results: 

Proposition 4.1. The negative polar cone of Cq defined by 

AAq := C° := { w G R^^ , (w, v) < Vv G Cq} , 
is equal to the proximal normal cone N(Q,q). More precisely, 

AAq = N(Q, q) = {- AijG,j(q), X,j > 0, Aj(q) > ^ = o} . 

Using the classical orthogonal decomposition with two mutually polar cones (see ^20j), the main 
equation (IIT]) becomes 

^ + PN(Q,q)(U(q))=U(q). (42) 

According to Proposition 3.3 in [2j, we know that for a Lipschitz map U, this differential equation 
is equivalent to the following differential inclusion: 

^+N(Q,q)9U(q). (43) 

For all {i,j), Dij is a convex function and belongs to C'^iUij), with 

:= {q G ]R2iV , |q»-qi|-r>0} 
and satisfies Assumption (|AOp with c = r/ Moreover it is obvious that 

Vq G , |Gij(q)| = V2 and \B^Dij iq)\ < -. 



As a consequence, Assumptions (lAip and ()A4p are satisfied. It remains to check the inverse 
triangle inequality (Assumption ()A5'P ). which is the aim of the following proposition. 

Proposition 4.2 (Inverse triangle inequality). 

There exists 7 > 1 such that for all q G Q, 



0!ij\Gij{q)\ <7 

(ij)G/{q) 



(ij)6^(q) 



where 

/(q) = {{i,j), i < j, Dij{q} = 0} and aij are nonnegative reals. 
Constant 7 can he fixed as follows 

N 

7 = 3\/2iV 



(-n(t). 



The next subsection is devoted to its proof. 
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Remark 4.3. Note the sign of coefficients aij. From a general point of view, this inequality is 
obviously wrong if these coefficients are just assumed real. Indeed, for N > 6, the cardinal of 
the set /(q) is strictly larger than 2N , which involves a relation between vectors Gjj(q). 

Remark 4.4. In \17^ , we have already proved such a result. The proof we propose here gives a 
smaller constant 7 and above all, it allows to better understand how the non-uniqueness of the 
Kuhn-Tucker multipliers appears. 



Proposition 12.91 Theorems 12.131 and 13.41 imply the next results. 

Proposition 4.5. The set Q C M^^, defined by is rj-prox-regular with a constant 



N 



6iV 



Theorem 4.6. Assume that U is Lipschitz. The Cauchy problem associated to '^^ well- 
posed in the set of the absolutely continuous functions and the related numerical scheme 12^) is 
convergent. 

Other numerical simulations based on this scheme, we refer the reader to [16 ^ \19 \ I17|. 
4.2.2 Proof of the inverse triangle inequality 

In order to prove Proposition 14.21 we are firstly going to show that the terms Xij appearing in 
the following equation : 

^A,,Gi,(q) =F 

i<j 

are bounded. 

Proposition 4.7. For all q G for all F G M^^, the following set 



A, 



JV(JV-l) 



XijGijiq) = F, Xij > 0, Xij = i/ Ai(q) > 



is uniformly bounded with respect to q. More precisely, 

VA G Aq^i?, Vi < j, Xij < |F| with a 



sin(f)' 



First we check that Proposition 14. 71 implies Proposition 14. 2[ Indeed, by Proposition 14. 7t we have 
for all {k, I) G /(q). 



N 



(ij)GJ(q) 



with 



sin(t)- 

Since |Gfci(q)| = v^, by summing coefficients Uki for all {k,l) G /(q) , we obtain 



«H|GH(q)| < V2 |/(q)| 

(fc,/)e/(q) 



(ij)G/(q) 
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In the monodisperse case, each disk has at most 6 neighbours. As a consequence, |/(q)| < 3N, 
and thus, 



J2 afcdGH(q)| < SV2N 
(fc,Oe/(q) 

It now remains to prove Proposition 14.71 



{*,i)eAq) 



Proof of Proposition \4- 7\ Suppose that set A^^p is not empty, we want to estimate the solutions 
A of the following system containing 2N equations: 

J2 A,,G,,(q) =F, Ai,- >0. (P) 

1 < j 
A,(q) =0 

By specifying the expression of vectors Gij(q), we can write the system concerning individual 

n 

i = i 

2 neighbour of io 



where 



T— _ / ^iio if J < ^0 and p • - e •• ici\ 



We assume here that the configuration q G Q is a cluster in the sense that the set IJi^i -^(li) ^) 
is connected by arcs. Otherwise, for the other configurations, it suffices to deal with the different 
clusters one by one, as the associated problems are independent from each other. We denote 
by ^ = {qi,..,qAr} the set of all positions. We are going to solve the problems (Pi) one by 
one. After having solved a problem (Piq), q^g is removed from A and problem {Pig) is deleted 
from the whole problem (P) , and the values of Ajjg are taken into account in terms Fj , for each 
person j neighbour of individual iq. The goal is to reduce yl to a singleton. 
The order of this algorithm is important. Less terms there are in problem (Pio), more easily 
solvable it is. First and foremost, we deal with people having only one neighbour (Case 1) 
because in that case the solution is trivial. When there is no single contact, we consider people 
with more neighbours. Since we want to have an upper bound for the terms Xjig, we will take 
care of controlling the angles between vectors ejpj^ and Ci^j^ if individuals ji and j2 are in contact 
with person zq (see Case 2). 

Case 1: elimination of single contacts 

Suppose that there exists qi G A such that person i has only one neighbour j. Problem (Pj) 
becomes 

Consequently, Xji = |Pj|. We remove q^ from A and (Pj) from (P). Then, we substitute Fj by 
Fj — XjiCij, which implies that the term \Fj\ is replaced with |Pj| + |Pj|. In this way, all single 
contacts are taken off, which allows us for example to completely deal with the case illustrated 
in Figure m If A is not reduced to a singleton after these eliminations, we consider the following 
case. 

Case 2: there are no single contacts 

We define C the convex hull of A. The set of its extremal points is denoted by E and the cardinal 
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Figure 4: Case of a chain. 
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Figure 5: Case 1. 

of E by p. We denote by P the boundary of C, which is a convex polygon whose vertices are the 
points of E. This polygon has p > 3 corners at least, otherwise we would rather consider Case 1. 
The angle sum equals to tt{p — 2). Therefore there exists an angle of P lower than 7r(l — |) < vr. 
Suppose that qj is the corresponding vertex. Because of the angle condition, individual i has two 
or three neighbours. Indeed, if individuals j and k are in contact with person i, the minimum 
angle between vectors Cij and Cik is equal to ^. 

Case 2 (a): individual i has 2 neighbours 

Using a rotation around q^ if needed, we can assume that the concerned configuration is repre- 
sented in Figure El In view of the choice of q^, 6* > y > ^ and Note that N >p>3. 
By writing the expression of Cji and e^j, we can solve Problem (Pj) after a simple computation: 

(Xji \ 1 ( ~ sin 6 cos 
/ sine 1^ 

Hence 

— , , — 1,1 

Aji < ^-^\Fi\ and < ^—\Fi\. 
sm V sm U 

We remove q^ from A and we replace \Fj\ (respectively with \Fj\ + V2\Fi\/ sin(27r/A^) (re- 
spectively \Fk\ + \Fi\/ sm{2TT / N)) because sin^ > sin(^). Problem (Pj) is also eliminated. 

Case 2 (b): individual i has 3 neighbours 

Using a rotation around q^ if needed, we can assume that the concerned configuration is rep- 
resented in Figure [H Because of the non-overlapping constraint, angles 9 and (f) satisfy the 
following double inequalities 

^<e<^. (44) 

AT - - 3 ^ ^ 

and 

2vr 27r , vr , , , 

<(h< . (45) 

AT 3 3 ^ ^ 

Note that N > p > 6. In fact, for N < 7 persons, we can see that the cluster A can be reduced 
to a singleton only by considering Case 1. 
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Figure 6: Case 2. 

In this case, the terms Xji, X^i and Xu are not unique (3 unknowns and 2 equations). We define 
a particular solution of Problem (Pi) (A°^, 0, A^J, where 
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sm( 



sine cose \ / Ff 
I ]\ Ff 



(46) 



To obtain all solutions, it remains to describe the kernel of System {Pi). We can easily show the 
next lemma. 

Lemma 4.8. The kernel of System {Pi) is generated by: 

sm{(p — 9) 
sin 6 
sin^ 

In addition, the signs of the coordinates can he specified 

sin(^— d) < — sin < 0' ^ ^ ^ smcf) < < 0- 



All the solutions of {Pi) take the form 







( \ 




( sm{4>-e) \ 









+ t 


sin 6 










y sin ^ J 



where t G 



Since we are only looking for the nonnegative solutions, the sign of sin 9 implies t > 0. Further- 
more, the signs of sin{(p — 9) and of sin (f) involve t < tmax where 



imax = min 
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sin(^ — 9)' — sin < 
Moreover, the following inequalities are satisfied: 



(47) 



— \/2,, — 1,1 — 2,, 

< and A^j < tmax sm ^ < ^ . 

sme* V3 



■' sin9 



So we remove q^ from A and we replace \ Fj\ (respectively and \Fi\) with |i<j|+\/2|-Fj|/ sin(27r/iV) 
(respectively + 2\Fi\/y/3 and + \Fi\/ sin(27r/A^)) . Problem (Pj) is also eliminated. 
Now we have removed one person from set A. If this set is not reduced to a singleton, we return 
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to Case 1 or Case 2. 



Conclusion 

By iterating this process, we solve Problem (P) with at most {N — 1) steps. We can conclude 
with the following lemma, whose the proof is a straightforward computation. 

3 

Lemma 4.9. At every step, the norm IFI is replaced at worst with — IFI. 

sin(f ) 

This lemma ends the proof of Proposition 14.71 □ 

Remark 4.10. A similar result can be proved in the polydisperse situation (where the radii 
are not assumed equal). Case 2 presents also more possibilities because the maximum number 
of neighbours that a person can have, depends on the radii ri. By denoting this number, a 
simple computation shows that 

IT 

riv < 



arcsm 



where rmin = minrj et rmax = maxr^. With these notations, it can be proved that in the 
polydisperse case (see \24^ ) 



2 /n 

VA G KqF, Vi < j, Ai,- < IFI b^ with b = — 

/ / TT \ / 27r 
mm sm , sm 



TT 



- 



Remark 4.11. The radii can be allowed to depend on time in a Lipschitz way. This new problem 
fits into the framework of sweeping process considered in this paper provided that 

r := inf inf ri(t) > 0. 
te[o,T] ie{l,...,N} 

Indeed Proposition ^-^ still holds with Ip{t, q), p < r by replacing 7 by another constant. 
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